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ABSTRACT 

Slow rotation and absence of strong mixing processes in atmospheres of chemically 
peculiar stars develop ideal conditions for the appearance of abundance anomalies 
through the mechanism of microscopic particle diffusion. This makes these objects 
look spectroscopically and photometrically different from their "normal" analogs. As 
a result, it is often difficult to accurately determine atmospheric parameters of these 
stars and special methods are needed for the consistent analysis of their atmospheres. 
The main aim of the present paper is to analyse atmospheric abundance and stratifica- 
tion of chemical elements in the atmosphere of the chemically peculiar star HD 103498. 
We find that two model atmospheres computed with individual and stratified abun- 
dances provide reasonable fit to observed spectroscopic and photometric indicators: 
Toff = 9300 K, \ogg = 3.5 and Tcff = 9500 K, \ogg = 3.6. It is shown that Mg has 
a large abundance gradient in the star's atmosphere with accumulation of Mg ions in 
the uppermost atmospheric layers, whereas Si demonstrates opposite behaviour with 
accumulation in deep layers. In addition, a detailed non-LTE analysis showed that 
none of Mg transitions under consideration is a subject of noticeable non-LTE effects. 
Comparing observed photometry transformed to physical units we estimated the ra- 
dius of HD 103498 to be between R = (4.56 ± O.77)i?0 for T^s = 9300 K, logy = 3.5 
and R = (4.39 ± O.75)i?0 for T^s = 9500 K, logy = 3.6 models respectively. We 
note that the lack of suitable observations in absolute units prevents us to uniquely 
determine the Toff of the star at the current stage of analysis. 

Key words: stars: chemically peculiar - stars: atmospheres - stars: individual: 
HD 103498 



1 INTRODUCTION 

As name indicates, Chemically Peculiar (CP) stars are un- 
usual in chemical composition in comparison to that of nor- 
mal stars with similar fundamental parameters. This family 
of stars resides on the upper main sequence and they con- 
tribute about 15% — 20% of early B to late F-type stars. The 
speciality of this group is that they hold the strong abun- 
dance anomalies in their atmospheres such as non-uniform 
horizontal and vertical distributions of chemical elements, 
and/or surface magnetic fields of different intensities. The 
type of peculiarity varies from star to star and it depends 
upon many factors like: stellar effective temperature, rota- 
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tional velocity, the presence of magnetic field, membership 
in binary stellar systems etc. Microscopic chemical diffusion, 
arising from a competition between radiative levitation and 
gravitational settling is believed to be the main process be - 
hind this phenomena, as firstly discussed bv lMichaudll|l970l ). 

HD 103498 is a characteristic example of an evolved 
star belonging to the co-called CP group of spectral type 
A (Ap) with strong overabundance of Si, Fe, and Cr in its 
at mosphere as follow s from the recent analysis undertaken 
bv lJoshi et all (|20ld ). The authors aimed on the pulsation 
analysis of the star's atmosphere and performed only classi- 
cal homogeneous abundance analysis based on scaled-solar 
metallicity model atmosphere. Their analysis did not pro- 
vide clear evidence whether HD 103498 belongs to roAp 
(rapidly-oscillating Ap) stellar group. They also noted that 
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the chemical properties and a high effective temperature 
of the star are distinguishable among known roAp stars. 
According to t he rec ent spectropolarimetric observations of 
lAuriere et all (1200711 this star has a weak longitudinal field 
< Bz > that varied between ±200 G with a period of 15.83 

days. 

In this work we extend the study of Ijoshi et al.l (|2010l ) 
and carry out complex and self-consistent analysis of the 
HD 103498 atmosphere. The main goal is to involve all the 
available observables to derive consistently fundamental at- 
mospheric parameters and abundance pattern of the star. 
We also look for stratification of chemical elements and 
attempt to construct a self-consistent atmospheric model 
which fits metallic line spectra, hydrogen line profiles, and 
energy distribution of the star. This approach where an it- 
erative process is used for the derivation of atmospheric pa- 
rameters was already successf ully applied in a number of re- 
cent investiga tions of normal (lRva,bchikova et al.ll2009|'l and 
pecu liar (e.g. IShulvak et all |2009| , I2OI0I : iKochukhov et all 
I2OO9I ) stars. The case of HD 103498 is of certain interest 
also because of its relatively high (as for roAp stars) tem- 
perature (Toff ~ 9500 K) at which such important spectro- 
scopic indicators as hydrogen lines become sensitive both 
to the atmospheric temperature and pressure stratification 
making it difficult to separate effects of T^s and logg. Thus 
other (photometric and/or spectrophotometric) observables 
must be simultaneously involved in the analysis to obtain 
precise results. 



2 OBSERVATIONS 

We used the average spectrum of HD 103498 which were ob- 
tained on February 2, 2009 (HJD 2454865.624-2454865.764) 
using the 2.56-m Nordic Optical Telescope (NOT) equipped 
with the Fibre-fed Echelle Spectrograph (FIES). The in- 
strument was configured to use the high-resolution mode 
at the resolving power of ii = 47 000 in the spectral region 
A A 3900-7270. The FIES spectra were acquired at the mag- 
netic phase range 0.37 ± 0.2. The abundance analysis of the 
program star is based on the average spectrum of time-series 
spectroscopy. The finer details about the observa t ions and 
data reduction procedure are given in IJoshi et all (|2010l l. 

We also m ade use of av ailable ph otometric observa- 
tions in Geneva (lRufeneij|l976l). Johnson (IHuchra &: WillneJ 
I1973I), Stromgren jHauck fc M crmiUiod' Il998l) 2M ASS 
i Skrutskie et al.ll2006l ). and A a (|Paunzen et alJl2005h sys- 
tems. 



3 ANALYSIS TOOLS AND METHODOLOGY 
3.1 Model atmospheres 

To perform the model atmospher e calculations we use d the 
recent version of the LLmodels (jShulvak et al.ll2004l ) stel- 
lar model atmosphere code. LLmodels is 1-D, hydrostatic, 
plain-parallel LTE code which accounts for the effects of in- 
dividual and stratified abundances. This program treats the 
bound-bound opacity by direct, line-by-line spectrum syn- 
thesis. The stratification of chemical elements is an input pa- 
rameter for the code and thus is not changing during model 
atmosphere calculation process. This allows to explore the 



changes in model structure due to stratification that was 
extracted directly from observations without modeling the 
processes that could be responsible for the observed inho- 
mogeneities. 

For every model atmosphere calculation the stellar at- 
mosphere is discretized into 120 layers between logrsooo = 
—8 and 2 with the higher point density in the region of steep 
abundance gradients to ensure accurate integration of radi- 
ation field properties and solution of other equations. The 
frequency-dependent quantities are calculated with variable 
wavelength step with a total of « 520 000 wavelength points 
spread over the spectral region where the star radiates most 
of its energy. For the co mputation of line a bsorptio n coef- 
ficien t VALD database jPiskunov et al. | [l995i : iKupka et al.l 
I1999I I was used as a source of atomic lines data. 



3.2 Classical abundance analysis 

The line identification in the observed spectrum of the star 
is based on the theoretical spectrum calculated for the entire 
observed spectral region using lines extr acted from VALD . 
The spectrum synthesis code SYNTH3 iKochukhovl (|2007l ') 
was used in all synthetic spectrum calculations. We com- 
pared synthetic and observed spectra to choose the least 
blended lines for the accurate abundance analysis using IDL 
package BinMagQ. 

The classical abundance analysis in the approximation 
of chemically-homogeneous atmosphere is based on the anal- 
ysis of equivalent widths of the lines, wh ich is performed us- 
ing the updat ed version (|Tsvmballll996l l of WIDTH9 code 
(|Kurucjll99l ). We adopted vsini = 12 km/s (obtained by 
the best fit to the observed unblended line profiles) and mi- 
croturbulent velocity ^ = 1.0 ± 0.2 km/s (deriv ed from nu- 
merou s lines of Cri, Crll, Fel and Fen) from [joshi et al.] 
|20l3)- As atmospheres of CP stars of type A are believed 
to be quiet in order to develop c hemical peculiarit ies, a non- 
zero microturbulence found in (jjoshi et al.ll2O10l ) probably 
results from the Zeeman broadening of spectral lines caused 
by a weak surface magnetic field. Since both in Joshi et al 
and in the present study we use non-magnetic spectra syn- 
thesis, ^ = 1.0 km/s was still needed to mimic the effect of 
the star's magnetic field and its value were kept constant 
during iterative abundance analysis described below. 



3.3 Stratified abundance analysis 

To study the stratification of chemical elements in the at- 
mosphere of the star we used the DDAFiT script which is 
an automatic procedure to fin d chemical abund ance gradi- 
ents from the observed spectra (|Kochukhovll2007l '). DDAFiT 
provides an optimization and visualization interface to the 
spectrum synthesis calculations. In this routine, the verti- 
cal abundance distribution of an element is described with 
the four parameters: chemical abundance in the upper at- 
mosphere, abundance in deep layers, the vertical position 
of abundance step and the width of the transition region 
where chemical abundance changes between the two values. 
All four parameters can be modified simultaneously with 
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the least-squares fitting procedure and based on observa- 
tions of unlimited number of spectral regions. This pro- 
cedure is succ essfully applied in a numb er of studies (e.g . 
Rvabchikova et al.., 2002. 2005, . ,200& : .Kochukhov et al.ll2009l : 



Calculation of first approximation 



Shulvak et al.ll2009l . boiij v 



3.4 Non-LTE analysis codes 

To analyse the non-LTE effects on the Mgl/Mgll ions we 
used the DETAIL and SU RFACE codes originally deve l- 
oped by Butler and Giddings ijButler ll984l : lGiddingsll98lh . 
Our calculations take into account the recent improvements 
in the atomic data for Mg. The extensive description of its 
model atom and non -LTE line formation are presented in 
IPrzvbilla etHI l|200ll ) and we refer the interested reader to 
the aforementioned paper. 



3.5 Self-consistent atmospheric modeling 

Using the analysis techniques outlined above, it is possible to 
reconstruct the element distribution profile of any chemical 
species for which the accurate atomic data exist. However it 
is to note that the empirical analysis of chemical elements 
stratification depends upon model atmosphere used, and the 
temperature-pressure structure of model atmosphere itself 
relies upon stratification. Consequently, the model atmo- 
sphere construction and stratified abundance analysis both 
are strongly coupled together. It is therefore necessary to 
use the iterative procedure which is described in different 
steps as follows: 

(i) selection of atmospheric parameters (first approxima- 
tion, i.e. initially estimated parameters must be close to 
those of the star to be analysed; they can be either taken 
from the literature or estimated using available photometric 
calibrations) ; 

(ii) analysis of individual spectral lines to determine the 
stratification of chemical elements; 

(iii) with the help of stratification found in previous 
step, new model must be calculated to provide improved 
temperature-pressure structure of the stellar atmosphere; 

(iv) comparison of the modeled energy distribution 
(or/and photometric colors) and hydrogen line profiles with 
the observed ones. The model input parameters (Toff, logjr) 
are then tuned to provide an agreement between theory and 
observations. 

(v) finally, the overall process from step 2 onwards must 
be repeated until stratification profiles of chemical elements 
and model parameters are converged. 

The mentioned above steps are schematically presented 
in Fig.[T] This procedure guarantees the consistency between 
the atmospheric model structure and abundances used for 
the calculation of synthetic line profiles. 

The final self-consistent, chemically-stratified model at- 
mosphere is expected to reproduce simultaneously the ob- 
served photometry, energy distribution, hydrogen line pro- 
files and metallic line spectra. Some re cent examples of 
apply ing t his procedure can be found in iKochukhov et akl 
l|2009l ) and lShulvaket al.l (|2009l . |201G| ). 



Determine stratification of chemicol elements 



Compute irrproved model atmosphere 



Comparison of modeled ond observed parameters 
(H — line profiles, ptiotometry and/or energy distribution, 
metollio line speotro) 



Check if stratification profiles and 
model parameters converged 



'not converged' 



'converged' 



Final results 



Figure 1. Steps used for the self-consistent models witli empirical 
abundance / stratification analysis . 



4 RESULTS AND DISCUSSION 

4.1 EYmdamental atmospheric parameters 

We used an iterative process in which the atmospheric pa- 
rameters Tcft , logg, and abundances are gradually improved 
by using different spectroscopic and photometric indicators. 
Using calibration of Stromgren indices and Geneva p hotom- 
etry as implemented in the package TEMPLoGG™'^ (|Kaiseij 
2006), as well as sim ultaneous fit to th e observed profiles 
of Ha and H/3 lines, IJoshi et al. 1 (2010') obtained the fol- 
lowing atmospheric parameters: Tcb — 9500 K, logg — 3.6 
employing ATLAS9 (|Kurucj|l992l . ll993l ') model atmosphere 
computed with scaled solar abundances [M/H] = -1-0.5 dex. 
We used this model as starting point in our analysis. Imple- 
menting the iterative procedure described in Sect. 13.51 two 
sets of parameters were obtained after five iterations as those 
providing, on average, a good agreement between observed 
and predicted spectroscopic and photometric observables: 
Tcff = 9300 K, logg = 3.5 and T^a = 9500 K, logg = 3.6. 
As will be discussed later, in some cases hotter models pro- 
vided better fit to certain photometric parameters, but the 
ionization balance and hydrogen lines were then subject of 
noticeable disagreement. 

The results of abundance analysis shows that the 
HD 103498 is Cr-F e rich star, and are consistent with those 
obtained before bv lJoshi et al] (|2010l ). Table [2] summarizes 
abundances obtained using two mod el atmosphe r es. Th e dif- 
ferences between abundances from Ijoshi et al.l (|20ld ) and 
ours for the T^g = 9500 K, logg = 3.6 model is likely due to 
the changes in atmospheric structure between solar-scaled 
abundance model and model atmosphere computed with re- 
alistic chemistry of the star. In addition, the number of lines 
used for abundance analysis is different between the two 
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Figure 2. Plots of individual abundances for Cr l/Cr ll(top) and Fe l/Fe II (bottom) lines as a function of excitation energy of the lower 
level for HD 103498 wih two different models; red filled circles - T^b = 9300 K, logg = 3.5; blue stars - T^s = 9500 K, logg = 3.6. 



studies which may also influence final results. Nevertheless, 
the agreement is very good within the estimated error bars. 
As expected, the model with cooler temperature statistically 
results on lower abundances as computed atomic lines be- 
come stronger with the decrease of Ten in this particular 
temperature regime. 

Figure [2] illustrates the ionization balance in the at- 
mosphere of HD 103498 as derived using lines of Fel/Fell 
and Crl/Crll. There is not much difference between the two 
models and the slopes are nearly the same, though model 
with Tcfi — 9500 K provides slightly better fit for first ions, 
while neutrals seem to prefer Toft — 9300 K. 



4.2 Stratification of Si and Mg 

Although HD 103498 is Cr-Fe rich star, these two elements 
do not demonstrate any signature of possible stratification. 
Abundances obtained from strong and weak lines of Fe and 
Cr were found to be nearly the same without noticeable 
patterns. Since lines with different excitation energies sta- 
tistically probe different atmospheric layers, the absence of 
peculiar abundance patterns among strong and weak lines 
clearly shows the absence of strong abundance gradients in 
the stellar atmosphere, observed in cooler CP stars (see, 
for example. Fig. 3 in Ryabchikova 2008). Nevertheless, we 
checked possible stratification with the DDAFiT package. 
We selected 17 lines of the neutral and singly ionised element 
for Cr and the same amount for Fe. The range of excitation 
energies were (0-3.45 eV) for Cri lines, (3.71-11.10 eV) for 
Crll lines, (0.86-3.39 eV) for Fel lines, and (2.90-10.68 eV) 
for Fell lines. The observed intensity varies from 11 to 300 
mA. Formally, for both elements we get a solution with small 
steep abundance gradient at log T5000 ~ — 1 which does not 
exceeds 0.5 dex, however, a difference between standard de- 



viation for homogeneous and stratified abundance distribu- 
tions is insignificant. For both elements, Cr and Fe, we get a 
small increase of abundance towards the upper atmosphere 
(except Cr for Tefi = 9300 K, logg — 3.5 ). Similar distribu- 
tions are observed in other CP stars of the same temperature 
range. Cr abundan ce gradient is very sma ll in HD 133792 
with Tcff = 9400 K (|Kochukhov et al.lbood') and practically 
absen t in HD 204411 with Tos = 8700 K (iRvabchikova et al.l 
2005h. Recent diffus ion calculations (|LeBlanc et al.l 20091 : 
Alecian &; stifti I2010I ) show that rather strong Cr and Fe 
abundance gradients with accumulation in deeper atmo- 
spheric layers should be observed in atmosphere of a star 
with TeS = 9300 — 9500 K having small magnetic field, 
and these gradients disappear at much higher temperatures. 
Strong horizontal magnetic field influence the diffusion pro- 
cess resulting in upward abundance gradient after abun- 
dance minimum at logrsooo ~ —1. In this case step- function 
approximation is not valid anymore, and its employment 
leads to smearing the possible gradients. However, no indi- 
cations exist in favour of the strong horizontal (toroidal?) 
component of magnetic field in HD 103498 as well as in 
HD 133792 and HD 204411. 

A detailed inspection of the spectra of HD 103498 
showed that only two elements are suitable for stratifica- 
tion analysis: Mg and Si. Both are represented by a sufficient 
number of strong and weak lines which are not blended much 
by nearby lines of others elements. These lines and their 
current parameters extracted from VALD are presented in 
Table [T] 

Using the DDAFiT package we derived stratification 
profiles of Mg and Si for two selected model atmospheres. 
Figure [3] demonstrates run of abundances as a function of 
atmospheric depth and Fig. [4] illustrates the resulting best 
fit to individual line profiles. We find magnesium to be 
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Table 1. A list of spectral lines used for the stratification analysis. 



Ion 


Wavelength 


i?i,(eV) 


log(9/) 


log7Stark 


Ref. 


Mgii 


4384.637 


9.996 


-0.790 


-4.07 


KP 


Mgii 


4390.514 


9.999 


-1.490 


-4.07 


KP 


Mgii 


4390.572 


9.999 


-0.530 


-4.07 


KP 


Mgii 


4433.988 


9.999 


-0.910 


-4.40 


KP 


Mgii 


4481.126 


8.864 


0.740 


appx 


KP 


Mgii 


4481.150 


8.864 


-0.560 


appx 


KP 


Mgii 


4481.325 


8.864 


0.590 


appx 


KP 


Mgi 


4702.991 


4.346 


-0.666 


-3.98 


LZ 


Mgii 


4739.593 


11.569 


-1.960 


appx 


KP 


Mgii 


4739.593 


11.569 


-0.660 


appx 


KP 


Mgii 


4739.709 


11.569 


-0.820 


appx 


KP 


Mgi 


5172.684 


2.712 


-0.402 


appx 


AZ 


Mgi 


5183.604 


2.717 


-0.180 


-7.99 


AZ 


Mgi 


5528.405 


4.346 


-0.620 


-4.46 


LZ 


Mgi 


5711.088 


4.346 


-1.833 


appx 


LZ 


Si II 


4130.872 


9.839 


-0.824 


-4.87 


BBCB 


Si II 


4190.707 


13.492 


-0.17 


-5.26 


MERL 


Si II 


5055.984 


10.074 


0.593 


-4.78 


SG 


Si II 


5056.317 


10.074 


-0.359 


-4.78 


SG 


Si II 


5669.563 


14.200 


0.28 


-5.53 


BBC 


Si II 


5688.817 


14.186 


0.13 


-5.50 


BBC 


Si II 


5868.444 


14.528 


0.400 


-5.36 


MERL 


Si II 


5957.559 


10.067 


-0.301 


-4.84 


SG 


Si I 


6155.134 


5.619 


-0.754 


-3.16 


K07 


Si II 


6347.109 


8.121 


0.297 


-5.04 


BBCB 



The columns give the ion identification, central wavelength, the 
excitation potential (in eV) of the lower level, oscillator strength 
(log{(;/)), the Stark damping constant ("appx" marks lines for 
which the classical approximation was used), and the reference 

for osci llator strength. 

AZ - lAnderson et al.l il96tl: BBC - iBlanco et all ll 19951): 
BBCB - iBerrv et al.l lll97ll): LZ - iLincke fc ZiegenbeinI 
lll97ll); KP - iKurucz fc Pevtrema rml l ll975D : K07 

|http://kurucz7har vard.edu/at oms/ 1400/ ME RL ^- 

IMatheron et al.l l l20oil'l : SG - ISmith fc Gallaghed l ll966f) . 
Stark damping constants are taken fro m VALD exc ept those 
marked by asterisk, which are taken from IWilke) 1I2OO3I') . 



strongly overabundant in upper atmosphere of HD 103498 
and under-abundant in deep layers. The abundance jump 
is very smooth, occupying a considerable region of atmo- 
spheric depths. It is to note that in case when all four pa- 
rameters of the assumed step-like abundance distribution 
were minimized simultaneously (see Sect. 13.31 for explana- 
tions), the solution was always with big error bars for the 
value of the abundance in deep layers. We thus also tried to 
derive stratification profile when abundances in these layers 
were fixed at constant values, from log(A'^/Aftotai) ~ —20 dex 
to log(Ai'/7Vtotai) ~ —10 dex (some examples are presented 
on Fig. [3)1. However, in all cases the fit to Mg lines did not 
change much as well as the common feature of strong abun- 
dance gradient between surface and photospheric layers stay 
still. It is important to stress here that any values of abun- 
dances below logrsooo ~ have no strict physical meaning 
because these atmospheric layers are not visible to us and 
thus cannot be modelled by comparing observed and pre- 
dicted line profiles. In other words, changing abundances in 



Table 2. LTE homogeneous atmospheric abundances of 
HD 103498. 



Ion This work Joshi et al. Sun 





t9300g3.5 


t9500g3.f 




Tl 


t9500g3.6 




Ci 


-4.20±0.15 


-4.11±0 


15 


2 


-4.11±0.23 


-3.65 





-3.81±0.22 


-3.79±0 


22 


2 


-3.92: 


-3.38 


Nai 


-4.93±0.16 


-4.82±0 


20 


2 


-5.09: 


-5.87 


Mgi 


-3.57±0.14 


-3.45±0 


14 


5 


-3.58±0.13 


-4.51 


Mgii 


-4.29±0.35 


-4.29±0 


35 


5 


-4.43±0.27 


-4.51 


Alii 


-6.21: 


-6.25: 




1 


-6.06: 


-5.67 


Si I 


-3.73±0.32 


-3.63±0 


32 


4 


-3.65±0.33 


-4.53 


Sin 


-3.94±0.36 


-3.99±0 


34 


7 


-3.64±0.47 


-4.53 


Can 


-5.87±0.09 


-5.79±0 


09 


2 


-5.91 


-5.73 


Ti II 


-6.46±0.17 


-6.35±0 


17 


25 


-6.45±0.16 


-7.14 


V II 


-8.17±0.18 


-8.07±0 


18 


2 


-8.14±0.18 


-8.04 


Cri 


-3.31±0.18 


-3.19±0 


19 


101 


-3.25±0.23 


-6.40 


Cm 


-3.36±0.22 


-3.32±0 


22 


230 


-3.31±0.24 


-6.40 


Mm 


-5.93±0.10 


-5.79±0 


10 


2 


-5.94±0.11 


-6.65 


Mnii 


-5.76±0.24 


-5.71±0 


24 


3 


-5.72±0.23 


-6.65 


Fei 


-3.04±0.16 


-2.92±0 


16 


105 


-2.98±0.20 


-4.59 


Fen 


-3.12±0.18 


-3.11±0 


18 


169 


-3.01±0.18 


-4.59 


Con 


-5.62: 


-5.57: 




1 


-5.58: 


-7.12 


Nil 


-5.33: 


-5.21: 




1 


-5.34: 


-5.81 


Ban 


-8.72±0.25 


-8.53±0 


25 


4 


-8.64: 


-9.87 


Ce II 


-9.16: 


-8.96: 




1 


-9.10: 


-10.46 


Prill 


-9.02±0.22 


-8.91±0 


22 


5 


-8.87±0.33 


-11.33 


Ndlll 


-8.51±0.16 


-8.40±0 


17 


12 


-8.41±0.16 


-10.59 


Smll 


-8.76: 


-8.58: 




1 


-8.70: 


-11.03 


Eull 


-8.94±0.14 


-8.72±0 


15 


2 


-8.85±0.12 


-11.53 


Gd II 


-8.63±0.19 


-8.47±0 


21 


2 


-8.72: 


-10.92 



Error estimates are based on the internal scattering from the num- 
ber of analysed lines, n. Abundance are given in log(Af/ A'total)- 
The last colum n gives the abundances of the solar atmosphere 
jAsplund et al.llioOSi') . 



Table 3. Non-LTE abundance corrections associated with each 
Mgi and Mgll transition used for stratification analysis. Col- 
umn 3 and 6 repres ents NLTE corrections for the star Vega 
llPrzvbilla et al.ll200ll'l . 



A log e A log e 

Mgi (hD103498) (Vcga) Mg II (hD10349s) (vega) 



4702.991 


-fO.Ol 


-1- 0.02 


4384.637 


- 0.00 




5172.684 


-0.03 


- 0.06 


4390.572 


- 0.00 


-1- 0.00 


5183.604 


-0.04 


- 0.13 


4433.988 


- 0.00 


+ 0.00 


5528.405 


-0.00 
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these layers has no effect on computed line profiles and hence 
the true abundances can not be recovered in principle. 

In addition we find that none of the Mg lines used 
for stratification analysis is a subject of strong deviations 
from the LTE which could affect final results. Since DE- 
TAIL/SURFACE package can only deal with homogeneous 
abu ndances , we u sed the original model atmosphere from 
J oshi et al.l (|2010l ) in order to estimate respective non-LTE 
abundance corrections A log e = log ejvLTE — log eltb , which 
are presented in Table [3] The largest non-LTE abundance 
correction was found for the Mgll 4481.126 A line and is 
-0.05 dex. 
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Stratification profile of Si was found to be similar to 
those frequently found in cooler CP stars, i.e. with accumu- 
lation in photospheric layers. Its shape did not change much 
between T^a = 9300 K and T^s = 9500 K models, only the 
abundance jump became more steep for the cooler one. 

The fact that Mg stratification is opposite to those of 
Si, i.e. that the Mg is pushed up by radiation much more ef- 
fectively is, of course, a very important result for the modern 
models of particle diffusion in stellar atmospheres, and the 
same behaviour of Mg was alr eady found in another pec uliar 
star HD 204411 (see Fig. 7 in lRvabchikova et aLllioOsI ). 

The obtained Mg distribution again can be explained 
qualitively by diffusion in the prese nce of very strong hor i- 
zontal magnetic field (see Fig. 3 in lAlecian fc StifdlioOTl ). 
In principle. Si gradient s are also predi c ted by diffu sion 
models, however, ne ither lAlecian fc Stif3 l|2007l . l2010h nor 
iLeBlanc etHI (|2009l) predict strong Si overabundance at the 
photospheric layers, that is commonly observed in all stars 
investigated up to now. In all diffusion calculation maximum 
of Si abundance is solar or slightly below the solar value. 

Formal distributions of Fe and Cr are also shown on 
Fig. |3l but because of relatively small abundance jumps 
(^ 0.5 dex) their distributions can be safely described as 
uniform (within the errors of standard deviations between 
observed and predicted line profiles). Therefore, even if 
present, such weak abundance gradients of Fe and Cr would 
in no way affect photometric and spectroscopic observables 
compared to the models computed with simple homogeneous 
distributions assumed in the present work. 



00 



-02 



-04 



z 



o 



-10 



-12 



'T3„=9500, log g = 3.5' 



-06 



-08 



\ 

V 

\ 



T^„=9500, log g = 3.6" 




-4 -2 


-6 -4 -2 




I09T5000 


Mg 


'free upper and lower abundance' 


Mg 


'fixed lower abundance = -20' 


Mg 


'fixed lower abundance - -10' 




'free upper and lower abundance ' 




"free upper and lower abundance " 


Cr 


"free upper and lower abundance " 



4.3 Energy distribution and hydrogen line profiles 

Unfortunately, the lack of observed spectrophotometric 
or low resolution spectroscopic energy distributions for 
HD 103498 did not allow us to uniquely derive its effective 
temperature. We therefore only attempted to use available 
photometric observations in different systems. FigureOillus- 
trates a difference between observed and computed photo- 
metric parameters for three atmospheric models. All photo- 
metric p arameter s were calculated using modified computer 
codes bv lKurucj (|l993l ), which take into account transmis- 
sion curves of individual photometric filters, mirror reflec- 
tivity and a photomultiplier response functi on. From the 
HIPP ARCOS parallax of vr = 3.37±0.56 mas (|van LeeuwenI 
I2OO7I ) the distance to the star is d = 297 pc, which results 
in redenning value of E(B — V) = 0.06 app lying extinction 
maps published bv lAmores fc Lepind (|2005h . In general, the 
two final models with Tefi = 9300 K and T^s = 9500 K are 
in reasonable agreement with observations except Stromgren 
ci and Geneva U — B for which we find the largest devia- 
tion of about 0.2 mag. As ci measures the amplitude of the 
Balmer jump, one would suppose that the true temperature 
or gravity of the star should be much higher than adopted 
in the present study. Indeed, increasing T^s up to 9700 K 
improves the fit to the ci significantly. However, in this case 
the fit to hydrogen lines become much worse. Last but not 
least, the temperature sensitive index b — y clearly requires 
lower temperatures. The same concerns gravity of the star: 
increasing log g also helps to get close to the observed value 
of ci, but then we find disagreement between observed and 
predicted hydrogen lines as well as in ionization balance be- 
tween Fe I /Fe 11 and Cr I /Cr 11. 



Figure 3. Stratification of Mg and Si in atmosphere of HD 103498 
derived using two models with Tgfj = 9300 K (left panel) and 
Teff = 9500 K (right panel) and different optimization settings 
for the step-like stratification profile. Full red line - optimization 
of all four parameters of the stratification profile of Mg; dashed 
blue line - lower abundance is fixed for Mg at log(Af/Aftotal) = 
—20 dex; dash-dotted green line - lower abundance is fixed for 
Mg at log(A''/A''total) = ~10 dex; dash-three dots magenta line - 
optimization of all four parameters of the stratification profile of 
Si; purple long-dashes — optimization of all four parameters of 
the stratification profile of Fe and finally dotted orange line - 
optimization of all four parameters of the stratification profile of 
Cr. 



To better understand the problem, we plotted ob- 
served photometric parameters transformed to absolute 
fluxes and compared them with model predictions on 
Fig. [6] To ob tained fluxes in physica l units, we used 
calibrations by iRufener fc Nicoleg (|l988l ) for the Geneva, 
Ivan der Bliek et al.l (Il996l ) for the 2MASS. and lBessell et all 
I 1998h for the Johnson photometry. As expected, infrared 
fluxes are not sensitive very much to the temperature and 
gravity of models, but the strong deviation in Johnson and 
Geneva ultraviolet U indices is clearly visible. 

The fit to the Balmer Ha and H/3 lines is presented on 
Fig. [7| Accounting for the stratification clearly helps to im- 
prove the fit to the observed line profiles compared to the 
spectra computed with only homogeneous abundances, but 
both spectra fail to fit the transition region between cores 
and wings of lines where theoretical profiles are systemat- 
ically too much wide. The same holds true for the cooler 
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0.6 - 



Wavelength 

Figure 4. A comparison between the observed and predicted 
profiles of Mg (top panel) and Si (bottom panel). Dotted black 
line — observations; dashed red line - synthetic spectra com- 
puted with homogeneous abundances (log(AfMg/JVtotal) = —3.87, 
log(A'^Si/Aftotal) = —4.95) (in case of Si for initial guess we used 
the homogeneous abundance which gives the best fit to line pro- 
files used for stratification analysis); solid blue line — synthetic 
spectra computed with the final stratified abundance distribu- 
tions shown by solid red lines in Fig. |3] 



model with T^s = 9300 K (not shown on the figure for the 
be tter view). Still, t he solar-scaled abundance model used 
bv lJoshi et al.l (|2010| ) seems to provide a best fit among all 
the models considered in this study, however one should not 
forget that this model is inconsistent with the true abun- 
dance pattern of the star: |joshi et al.l l|2010l ) used ATLAS9 
model with metallicity [M/H] = -1-0.5 dex, but the true iron 
overabundance of the star is « -1-1.5 dex compared to the 
Sun. Obviously, using now respective homogeneous model 
with [M/H] = +1.5 dex would immediately result in wider 
profiles of hydrogen lines (due to increase in number of free 
electrons and thus Stark broadening) no longer providing a 
fit to line profiles. 

From the analysis presented above its clear that there 
is a fundamental disagreement between temperatures of the 
star suggested by some photometric indicators and hydrogen 
line profiles. Fitting energy distribution require systemati- 
cally hotter temperature than those followed from fitting 
the Ha and H/3 lines. We thus speculate that there could 
be few reasons for this that should be checked in future in- 
vestigations. First of all, except Mg and Si other elements 
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Figure 5. Difference between observed and predicted colors of 
different photometric systems. Theoretical colors were reddened 
applying E{B - V) = 0.06. 

can be stratified in the atmosphere of HD 103498. If so and 
(as usually found) the stratification profiles are similar to 
that of Si, then the concentrations of free electrons in upper 
atmospheric layers could becomes smaller than the one of a 
present study, thus resulting in narrower profiles of hydro- 
gen lines. Unfortunately, as stated above, no good sets of 
lines for accurate stratification analysis could be found for 
elements over than Mg and Si in the spectra of HD 103498. 
Also, the normalization of hydrogen lines in echelle spectra 
is a complicated task and some systematic uncertainties may 
still influence fitting results. 

Finally, comparison of observed photometry calibrated 
in absolute units and model fluxes allowed us to estimate 
the radius of the star. With parallax tt = 3.37 ±0.56 mas we 
obtained R = (4.56 ± O.77)i?0 for Tcfi = 9300 K, logg = 3.5 
model and R = (4.39±0.75)i?Q for Tcft = 9500 K, logg = 3.6 
model respectively. 



5 CONCLUSIONS 

In this work we carried out detailed atmospheric abundance 
and stratification analysis of the chemically peculiar AO star 
HD 103498 , whic h is an extension of analysis presented in 
iJoshi et al.ll|2010l ). Using available spectroscopic and photo- 
metric parameters, as well as up-to-date model atmosphere 
and spectra synthesis codes we attempted to construct a self- 
consistent model atmosphere of the star and to determine 
fundamental atmospheric parameters iteratively by compar- 
ing various observables with the model predictions. Our 
main results are summarized in following: 

• We find that, generally, two model atmospheres com- 
puted with individual and stratified abundances provide rea- 
sonable fit simultaneously to spectroscopic and photometric 
indicators: T^h = 9300 K, log 3 = 3.5 and T^a = 9500 K, 
logg = 3.6. The later model ha s the same parameters as the 
one used in IJoshi et al.l (|20ld ). but, in contrast, includes 
realistic chemistry of the star. 

• The abundance analysis demonstrate strong overabun- 
dance of Fe and Cr relative to the Sun, and o ur results are 
in excellent agreement with iJoshi et al.l(|201Gl ). 

• Using simplified step-function approximation of strati- 
fication profile we derived depth-dependent distributions of 




Figure 6. Comparison between observed photometric parameters calibrated in absolute units and theoretical fluxes convolved with 
FWHM=5 A Gaussian. Theoretical fluxes were reddened applying E{B — V) = 0.06. 
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Figure 7. Comparison between the observed H/3 (left panel), Ha (right panel) observed and s ynthetic line profi les. Solid black line 
- observed profile; dashed red line — synthetic spectra cumputed with homogeneous abundances jjoshi et al.ll20l"ch : dash-dotted green 
line - synthetic spectra cumputed with homogeneous abundances (this work); dash-three dots blue line - synthetic spectra computed 
with stratified abundances. In both cases the final adopted model corresponds to Tcfi = 9500 K, logg = 3.6 

. The inset shows the zoomed part of the profiles around the core of line. 



Mg and Si. We show that Mg has a large abundance gradi- 
ent in the atmosphere of HD 103498, with accumulation of 
Mg ions in the uppermost atmospheric layers. Distribution 
of Si demonstrates opposite behaviour with accumulation of 
Si ions in deep atmospheric layers. This indicates that Mg 
is pushed up by radiative forces much more efficiently than 
Si. These empirical findings are of certain interest for the 
modern theory of particle diffusion in stellar atmospheres. 

• Although the star is Cr-Fe rich, these two elements do 



not show any signature of noticeable stratification in its at- 
mosphere. 

• We find that none of Mg transitions used for stratifi- 
cation analysis is a subject of noticeable non-LTE effects. 
The largest non-LTE abundance correction was found for 
the Mg 114481.126 A line and is -0.05 dex. 

• Comparing observed photometry in Stromgren, 
Geneva, 2MASS, and Johnson filters transformed to phys- 
ical units with model fiuxes we estimated the radius of 
HD 103498 to be i? = 4.56 ± 0.77 Rq for T^a = 9300 K, 
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logs = 3.5 model and R = 4.39 ±0.75 Rq for T^g = 9500 K, 
logg = 3.6 model respectively. 

• Comparison of observed photometric parameters and 
hydrogen line profiles for different model atmospheres 
showed that such indices as ci and U — B systematically 
require hotter models and Ha, H/3, and b — y index cooler 
models relative to the initially assumed temperature Tes = 
9500 K. Unfortunately, the lack of suitable spectrophotomet- 
ric or low resolution spectroscopic observations in absolute 
units forbids us to uniquely determine the T^ff of the star. 
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